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We present improved measurements of the branching fractions and CP-asymmetries in the B" tt+tt , 
— > TT^Tr", and B+ — > p"'"p'' decays, which impact the determination of a. We find 

= -0.68 ±0.10 ±0.03 
C+- = -0.25 ±0.08 ±0.02 
0°^ = -0.43 ±0.26 ±0.05 
23(5" tt'VO) = (1.83 ± 0.21 ± 0.13) X 10"^ 

for B — ► 1T1T decays, and 

B{B+-*p+p°) = (23.7 ± 1.4 ± 1.4) X 10-*^ 
fhip'^P^) = 0.950 ± 0.015 ± 0.006 
app = {92.4l6;0)o 

for B ^ pp decays. 

The combined branching fractions of iJ — > ii'i(1270)7r and B i<'i(1400)7r decays are measured for the first 
time and allow a novel determination of a in the — > ai{1260)^7rT decay channel. We obtain 

^ ii:i(1270)+7r- + E'i(1400)+7r-) = {3.ll°;^) X lO''^ 

BiB+ /ri(1270)''7r+ +ft:i(1400)°7r+) = (2.81?;?) X lO''^ 

aai,r = {79±7±11)°. 

These measurements are performed using the final dataset collected by the BaBar detector at the PEP-II 
B-factory. 



1. Introduction 

The primary goal of the experiments based at the 
B factories is to test the Cabibbo-Kobayashi-Maskawa 
(CKM) picture of CP violation in the standard model 
of electroweak interactions [l| . This can be achieved 
by measuring the angles and sides of the Unitarity 
Triangle in a redundant way. 

An effective value aeff for the CKM phase a = 
arg(— V^j;V^^/y„^V^*j) can be extracted from the time- 
dependent analysis of B meson decays dominated by 
tree-level b —> uud amplitudes, such as B'^ — > tt'^tt", 
B° p+p-, B" p^TT^, and ai(1260)±7rT. 
The current average values of a, a = (92 ± 7)° 
and a = (89j;42)° [S], obtained with different statis- 
tical techniques, are based solely on the analysis of 
B — > TTTT, B pp, and B ^ pn decays. 

The measurement of the angle a has witnessed 
significant progress over the past year. The follow- 
ing sections are organized as follows: a brief intro- 
duction on the experimental technique is given in 
Sec. [H the summer 2008 update of the measurement 
of the time-dependent CP-violating asymmetries in 
TT+vr^ decays and of the branching fractions 
(BFs) of i?° TT^TT^ decays is reported in Sec. [31 
Sec. m describes the 2009 update of the BF measure- 
ment of B^ — > decays, and its impact on the 
precision of the determination of a; in Sec. [5l we in- 
troduce the first measurement of B — + /'iri(1270)7r and 
B _ft'i(1400)7r decays and a new determination of a 



in B° ai(1260)±7r=F decays. 



2. Experimental techniques 

The interference between the direct tree decay 
(which carries the weak phase 7) and decay after 
mixing (which carries a weak phase 2/3) results 
in a time-dependent decay-rate asymmetry that is sen- 
sitive to the angle 2/3 + 27 = 27r — 2a. 

At the asymmetric-energy e+e~ B-factory PEP-II, 
running at a center of momentum (CM) energy of 
10.58 GeV, a BB pair is coherently produced in the 
decay of a T{AS) resonance. The resulting BB system 
has a boost f3j ~ 0.56 with respect to the laboratory 
frame. By means of this experimental device it is pos- 
sible to measure the decay vertex displacement Az 
of the two B mesons in the event, and hence their 
proper-time difference Atmeas ~ f^- 

One of the B mesons (Brec) is fully reconstructed 
according to the final state of interest. In order to 
study the time-dependence of the decay rates, it is 
necessary to measure the proper-time difference At 
between the two B mesons in the event and to identify 
the flavor of the other B-meson (Btag)- The flavor 
and the decay vertex position of Btag sue therefore 
identified from its decay products. 

The decay-rate distribution for B'^ (P°) decays to 



2 
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a CP-eigenstatc, such as tt+tt , is given by: 



dN 
dAt 



HAt|/r 
17 



1 - qtag 



C coa{AmdAt) 



Ssm{AmdAt) 



(1) 



where r = (1.536±0.014) ps 3 is the mean B hfetime, 
Arud = (0.502 ± 0.007) ps"^ is the B° - 5° mixing 
frequency and qtag = +1 (—1) if the Btag decays as 
a (B'^)- The parameters S* and C describe mixing- 
induced and direct CP- violation, respectively, and are 
defined as: 



S 



2ImA 



c 



l + |A|^ 



with A = where q/p is related to the B^ — 5° 
mixing, and A {A) is the amplitude of the decay of a 

B° (B*^) to the final state under study. If only the tree 
amplitude contributes to the decay, S — sin(2a) and 
C = 0. However, b — > uud transitions receive sizeable 
contributions from penguin (loop) amplitudes, which 
carry different strong and weak phases. This contribu- 
tion can result in non-zero direct CP- violation C 7^ 
and modifies S into 



S' = sin(2aeff)Vl-C2. 



(2) 



The angle acs coincides with a in the limit of van- 
ishing penguin contributions. In order to constrain 
Aa = a — ttoff , techniques based on the SU(2) isospin 
symmetry (for decays to a CP-eigenstate, such as 
P" 7r+7r~, /9+yO~) or the SU(3) approximate flavor 
symmetry (for decays to a non CP-eigenstate, such 
as P° — > p^TT^, ai(1260)^7r^) have been devised, and 
are discussed in the remaining of this paper. 

A neural network based tagging algorithm Q is 
used to determine whether the Btaa is a P or a 
Events are separated according to the particle content 
of the Btag final state into events where there are lep- 
tons, kaons and pions, for a total of seven mutually 
exclusive categories. The performance of the tagging 
algorithm is characterized by the efficiency etag in the 
determination of the flavor of Btag and by the mistag 
probability uj, and depends on the tagging category. 
The At distribution of Eq. [1] is convolved with a de- 
tector resolution function, which differs for signal and 
background, and is parameterized as a triple Gaus- 
sian. Dilution from incorrect assignment of the flavor 
of Btag is also taken into account: 



dN 



\At\/r 



dAtr 



4r 



gtQg(l - 2uj) 



5P(Ai„ 



1 - qtagAu! - 

Ccos(AmdAt) 

S am{AmdAt) 
-At), 



where (1 — 2uj) is the dilution factor, Auj is the dif- 
ference in mistag probabilities Aw = ajgo — togo and 
R is the resolution function. The parameters of the 
resolution function are obtained from a fit to a large 
sample of fully reconstructed P decays, as in [S'j, and 
are free to differ between tagging categories. 

The analyses of the two-body and quasi-two-body 
decays described in the remaining of this paper rely 
on a common strategy for the suppression of the 
continuum e+e^ qq background {q ~ u,d,s,c), 
which represents the most abundant source of back- 
ground. Two kinematic variables, the energy substi- 
tuted mass ruEs — — and the energy differ- 
ence AP = Eb — Vs/2, where -y/s is the e+e" CM 
energy and the four- momentum (Eb,Pb) of the B 
meson is defined in the CM frame, allow to discrimi- 
nate correctly reconstructed P candidates (for which 
the distribution of mss peaks at the P-meson mass 
and that of AP peaks at zero) and fake candidates 
resulting from random combination of particles (for 
which niES follows a phase-space distribution and AP 
is approximately fiat). Topological variables provide 
further distinction between the jet-like shape of con- 
tinuum events and the more isotropic P decays, and 
can be combined into multivariate classifiers, such as 
neural network and Fisher discriminant, to enhance 
the discriminating power. The signal and background 
yields and CP asymmetries are extracted via an ex- 
tended unbinned maximum-likelihood (ML) fit to the 
data. 



3. Isospin analysis of B — s- tttt decays 



3.1. 



TT+TT" 



In the TTTT system the penguin pollution is greatest. 
The tree (P) and penguin (P) amplitudes each con- 
tribute, with different weak {(f)) and strong (S) phases, 
with comparable magnitude. Direct CP violation, 
which is given by ^cp = 2sin0sin(5/(|T/P| + |P/T| + 
2cos(/)COS(5), can, therefore, be within observational 
reach. 



P° TT+TT^ decays are analyzed with the full 
BaBar dataset of 467 ± 5 milhon PP pairs 0]. A 
simultaneous ML fit to the tt+tt^, n^K~, K^tt~, 
and K'^K~ final states is performed. K — w sepa- 
ration is obtained by particle-identification (PID) ob- 
servables (the Cherenkov angle Qc in the DIRC 
and ionization-energy loss dE/dx in the tracking de- 
(|3^ vices ^]). Additional separation between the final 
states under study is achieved from AP: since the P 
meson is reconstructed from two oppositely charged 
tracks that are both given the pion mass hypothe- 
sis, each charged K in the final state results in a 
AP displacement of about —45 MeV. We extract 
1394 ± 54 signal events. From the time distribution 
of P*^ TT+TT^ decays a non-zero mixing-induced 
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Figure 1: 5"^- and C+- in B° 



the central val- 



ues, errors, and confidence-level (CL) contours, calculated 
from the square root of the change in the value of —2 In £ 
compared with its value at the minimum [0] . The system- 
atic errors are included. The measured value is 6.7a from 



the point of no CP violation (St-w ~ and 



0). 



CP violation asymmetry S+- = -0.68 ± 0.10 ± 0.03 
is observed vifith significance 6.3cr as shown in 
Fig. [TJ A non-zero direct CP violation asymmetry 
C+- = -0.25 ± 0.08 ± 0.02 is also extracted with sig- 
nificance 3. Oct [i]. 



3.2. B° 



77 TT 



The -B'^ — > tt'^tt" decay is formed from pairs of 
7^0 candidates, where one of the photons can 

eventually be reconstructed from two tracks coming 
from a photon conversion 7 e+e^ inside the elec- 
tromagnetic calorimeter. 

The yield and the flavor tag- and time-integrated 
CP asymmetry A'fjp = — are obtained from a 
ML flt to the kinematic variables AE and tues and 
the output of a neural network NN computed from 
event-shape variables, as well as the output of the B- 
flavor tagging algorithm. The background model ac- 
counts for correlations between NN and niEs- We ob- 
serve 247±29 signal events (corresponding to B{B''^ 
7r°7r°) = (1.83 ± 0.21 ± 0.13) x 10"^) and extracts 
C°° = -0.43 ±0.26 ±0.05 3- Since no reliable vertex 
information is extracted, 5™ can not be determined. 



3.3. Isospin analysis ot B ^ nn decays 

The rates and CP asymmetries of tt+tt^ and 

tt'^tt^ decays are combined with the results for 



the TT^TT^ mode in a model-independent isospin 

analysis Q. Under the isospin symmetry, B nn 
amplitudes can be decomposed in isospin / = (Aq) 
and 1 = 2 {A2) amplitudes. By virtue of Bose statis- 
tics, 1=1 contributions are forbidden. The following 
relations hold 



(l/^/2)A^ 



A2 - Ao, 



A"° = 2A2 + Ao, A+" = 3A2 



(4) 



(5) 



where A'l {A ' ) are the amplitudes of B [B) decays to 
the tt'tt-' final state. This yields the complex triangle 
relations: 



1 



V2 



A+- = A+" - 



^A^-=A-'-^\ 
V2 



(6) 



(7) 



Tree amplitudes receive contributions from both Aq 
and A2, while giuonic penguin diagrams are pure / = 
amplitudes and do not contribute to B^ tt+tt^ 
amplitudes. Possible contributions from electroweak 
penguins (EWP), which do not obey SU(2) isospin 
symmetry, are assumed to be negligible and are there- 
fore ignored. Under this assumption, |^^°| = \A \ 
(a sizeable contribution from EWPs would result in 

^ \A ^\ and would be signalled by an evidence 
of direct CP violation in B+ tt+tt" decays). If 

j4+° and A are aligned with a suitable choice of 
phases, the relations © and ([7]) can be represented in 
the complex plane by two triangles (Fig. ^ , and the 

phase difference between A'^ and A'^ is 2Aa. 




Figure 2: Triangles in the complex plane describing the 
isospin relations Eq. ^ and Eq. ([7]). 

Constraints on the CKM angle a and on the pen- 
guin contribution Aa are obtained from a confidence 
level (CL) scan over the parameters of interest, a 
and |Aa|. Assuming the isospin-triangle relations ([6]) 



and (O and the expression ^ 
plitudes (A+°,A 



,A°°A'^~ 



a X for the five am- 
A ) is calculated from 
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Table I Summary of the input to the isospin analysis of 
the TTTT system 6', lOj. 



Mode 


B(xlO-^) 


C 




5.5 ±0.4 ±0.3 
5.02 ±0.46 ±0.29 
1.83 ±0.21 ±0.13 


-0.25 ±0.08 ±0.02 
(-0.03 ±0.08 ±0.01) 
-0.43 ±0.26 ±0.05 



the measurements summarized in Table HI and mini- 
mized with respect to the parameters that don't enter 
the scan. The 1 — CL values are then calculated from 
the probability of the minimized . 

The results of the isospin analysis are shown in 
Fig. [3] and Fig. [H Aa is extracted with a four-fold 
ambiguity, which can be graphically represented as 
a flip of either triangle around A^''\ An additional 
two-fold ambiguity arises from the trigonometric re- 
lation = sin(2aoff) VI — • This results in 
a global eight-fold ambiguity in the range [0, 180]° on 
the extraction of a. A value Aa < 43° at 90% CL is 
obtained, which dominates the uncertainty on a @. 
Considering only the solution consistent with the re- 
sults of global CKM fits, a is in the range [71, 109]° 
at the 68% CL :6]. 

The limiting factor in the extraction of Aa is the 
knowledge of \A^^\ and \A |, which is severely lim- 
ited by the available statistics. A significant increase 
in statistics is therefore required in order to perform 
a precision measurement of a in this channel. A mea- 
surement of S*^^ , which would aid resolving some am- 
biguities on a, can only be addressed with Super B 
factory luminosities [ll| . 




la - a J (degrees) 



Figure 3: Projection of the 1 — CL scan on Aa for the tttt 
system 0|. 



4. Isospin analysis ot B ^ pp decays 

With respect to B ^ tttt decays, B ^ pp decays 
have a more favourable penguin to tree amplitude 
ratio. Moreover, the BF for B^ — > p'^p~ decays is 
greater than that for 5° tt+tt" decays by a fac- 
tor of ~ 5 Finally, the B^ — > decay can 
be reconstructed from a final state consisting of all 
charged tracks, with enough efficiency to allow for a 
measurement of S"™ with the present statistics [Tsl ]. 
Despite these many advantages with respect to the 
isospin analysis of tttt decays, the pp system exhibits 
some potential complications. 

In B° P~^P~ transitions, a pseudo-scalar par- 
ticle decays into two vector mesons. Three helicity 
states {H = 0, ±1), with different CP transformation 
properties, can therefore contribute to the decay [l3 |. 
The H = state corresponds to longitudinal polar- 
ization and is CP-even, while the transverse polar- 
ization states H — +1 and H — —1 (which are su- 
perpositions of S-, P-, and D-wave amplitudes) have 
not a definite CP-eigenvalue . Isospin relations similar 
to Eq. ^ and ([7]) hold separately for each polariza- 
tion state. The analysis of the angular distribution of 
B^ p~^ p~ decays allows to determine the longitudi- 
nal polarization fraction f^: 

1 d^r 

-Fl H- 4/l cos^ 6li cos^ 6*2 

1 c(cost^iacost^2 

±(1-/L)sin2 0isin2 02, (8) 

where 9i (6*2) is the angle between the daughter tt" 
and the direction opposite to the B direction in the 
p+ {p~) rest frame, as shown in Fig. [S] Since ex- 
perimental measurements have shown the decay to be 
dominated by the longitudinal, CP-even polarization, 
it is not necessary to separate the definite-CP contri- 
butions of the transverse polarization by means of a 
full angular analysis. 

A second complication arises because the p mesons 
have finite width, thus allowing for the two p mesons 




""O 50 100 150 

a (degrees) 

Figure 4: Projection of the 1 — CL scan on a for the tttv 
system 0. 
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Figure 5: Definition of the 
decays [l^ . 



^1 and 62 angles in B 



PP 



in the decay to have different masses. Since the Bose- 
Einstein symmetry does not hold, the wave function 
of the pp system can be anti-symmetric, and isospin 
1=1 amphtudes are allowed, breaking the isospin 
relations Eq. (O and ^ The stability of the fit- 
ted CP- violation parameters against the restriction of 
the TTTT invariant mass window used to select the p can- 
didates shows however that possible isospin violation 
effects are below the current sensitivity. 



4.1. B+ 



super- 
-1 



The decay analysis has been updated 

using the final BaBar dataset of 424 fb~^ [16| 
seding the previous analysis based on 211 fb~^ [1 
An analysis of the angular distributions of _B+ - 
decays is performed. The signal yield and longitudinal 
polarization fraction is extracted via a ML fit to the 
kinematic quantities m^S: Ai?, the output of a neu- 
ral network NN based on event-shape variables, the 
mass of the p+ and p^ candidates, and the cosines of 
the helicity angles 9p+ and 9pO , where 6p+ [OpO ) is the 
angle between the daughter tt*' (tt^) and the direction 
opposite to the B direction in the /9+ [p^) rest frame. 

Improvements have been introduced on the charged 
particle reconstruction and on the background model, 
which takes into account correlations between NN , 
the cosine of the helicity angle, and the tttt invariant 
mass for each p meson in the final state. The mea- 
sured BE increases from (18.2 ± 3.0) x 10"^ [H up 
to (23.7± 1.4± 1.4) X 10-6 The longitudinal po- 
larization fraction is /l = 0.950 ± 0.015 ± 0.006 
The measured direct CP- violation asymmetry Ac p = 

r(B-^p-pO)+r(s+-p+p") ^cp - ± U.U55 ± 

0.010, which is consistent with 0. This result indi- 
cates that the contribution from EWPs is negligible, 
and the i sosp in analysis holds within an uncertainty 

of 1 - 2° m. 




-20 



20 
Aa^-agff(deg) 



Figure 6: Projection of the 1 — CL scan on Aa for the pp 
system [l6| . 



^ 1 



0.5! 







ii (1 

I \ CL=68% \ 1 
\ \CL=90%/ / 


I I 







50 100 150 

a (deg) 



Figure 7: Projection of the 1 — CL scan on a for the pp 
system p^ . 



The BEs, longitudinal polarization fractions, di- 
rect and mixing-induced CP violation asymmetries 
for B ^ pp decays are used as input to the isospin 
analysis, and are summarized in Table [III The BE's 
of P^ P^ P^ and P° — *■ p^ p~ are now very sim- 
ilar and much higher than that for the B'-^ p^ p^ 
penguin transition. As a consequence, the isospin tri- 
angles do not close, i.e. \A+-\/^+ < 
This results in a degeneracy of the eight-fold ambi- 
guity on a into a four-fold ambiguity, corresponding 
to peaks in the vicinity of 0°, 90° (two degenerate 
peaks), 180°, as shown in Eig. [Hland Eig. [T] A value 
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Table 11 Summary of the input to the isospin analysis of the pp system [l^ . [T^ . [l^ . 



Mode 




fL 


c 


S 


p+p- 
p+p" 


25.5±2.li|;5 
23.7 ± 1.4 ± 1.4 
0.92 ±0.32 ±0.14 


0.992 ± 0.024l° E;26 
0.950 ±0.042 ±0.006 
0.75lg:ll ± 0.04 


0.01 ±0.15 ±0.06 
(0.054 ± 0.055 ±0.010) 
0.2 ±0.8 ±0.3 


-0.17 ±0.201^6 
0.3 ±0.7 ±0.2 



-1.8° < Aa < 6.7° at 68% CL is obtained. Consid- 
ering only the solution consistent with the results of 

1-6.0 
-6.5- 



global CKM fits, a = 92.4+^;|?. The precision on a is 



now at the level of 5%. 



5. ai(1260)±7rT 

It is possible to extract a from B decays to final 
states that are not CP-eigenstates [l^, such as B'^ 
ai(1260)^7rT decays. The relevant amplitudes are: 

A+ = A{B^ a+Tj-), A+ = ^ ar^+)> (9) 

yl_ EE ^ ar7r+), 3_ = A{B^ a+7r-).(10) 



The time distribution for this decay mode is given by: 

e-|A*l/r 



dAtr, 



(1±A 



CP) 



At 



1 - qtag^^ 



iRiAtr, 



iS±AS) sin(AmdAi) ■ 

(C ± AC) cos(A77idAi) 
-At), 



where Acp is the time- and flavor- integrated charge 
asymmetry, and 



S-i AS* 



C± AC 



2Im (e-2*^y4=pAi) 



(11) 

(12) 



The measured CP- violation parameters for P" - 
ai(1260)^7r^ decays are summarized in Table lllll [2C 



Table III Values of the CP- violation parameters used as 
input to the calculation of the bounds on |Aa| [20| . 



Acp 


-0.07 ±0.07 ±0.02 


S 


0.37 ±0.21 ±0.07 


AS 


-0.14 ±0.21 ±0.06 


C 


-0.10 ±0.15 ±0.09 


AC 


0.26 ±0.15 ±0.07 



In analogy to the tt+tt case, where 
2acff = arg [e''^"^ A{B° tt+tt-)A*{B° Tr+Tr")] 
it is possible to define two quantities and a^: 



± ^arg [e-2''^A±Ay , 



(13) 



which are related by the phase 5 = arg[A-^A*_] to the 
measurable quantities: 

arg [e-^'f^A^A*^] (14) 
St AS 



— arcsm 



VI-(CtAC)^ 



(15) 



In the limit of zero penguin amplitudes, S coincides 
with the strong phase difference between the tree 
amplitudes contributing to P° ai(1260)"'"7r~ and 
P" ai(1260)^7r+ decays. An effective value a^s 
for the weak phase a is then obtained as the aver- 
age Ooff = ^ {o^ts ~^ '^cs) "^ith an eight-fold ambigu- 
ity 

It is possible to apply arguments based on the 
approximate SU(3) flavor symmetry to set bounds 
on |Aq;|. The following ratios of CP-averaged rates 
of AS = and AS* = 1 transitions are calcu- 
lated, that involve the same SU(3) flavor multiplet 
as ai(1260) [U, such as P° ^ ai(1260)-/-i:+, P" ^ 



K 



lA 



TT , P+ ai(1260)+ii:°, and P 

P, 
P° 

Pi 



r+- 



— T ^ ■ 



f}-BiB"^a-^+) ' 

X^f^BiB + ^a+K") 
flBiB°^a-^+) ■ 



(16) 
(17) 
(18) 
(19) 



The bounds are effective because the penguin con- 
tribution is CKM enhanced by 1/A = | | / 1 Krf | in 
AS" = 1 decays with respect to AS* = modes. The 
following inequalities involving (a^ — a) hold: 



cos2(Q!^fj -a)> 
cos2(a^ — a) > 



1 - 2P5. 



I- A 



±2 



2P+ 



CP 

f 
± 



1-A 



±2 
CP 



(20) 
(21) 
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where A^p are the direct CP asymmetries 



A 



\A±\ 



CP — i-r 



(22) 



The above relations set a constraint on (a^f 
Bounds on |Aa| are then derived from |Aa| < 
^cff-"l)/2- 



-a). 

4s- 

The BFs oi B ^ ai(1260)7r and B aiil260}K 
decays have been measured in the last few years [20[. 
The measurement of the missing piece of input, the 
BFs of B ^ ifi(1270)7r and B Js:i(1400)7r decays, 
is described in the following section. 

5.1. B Ki(1270)7r, B Ki(1400)7r 

The KiA meson (the SU(3) partner of the ai(1260) 
meson) is a nearly equal superposition of the phys- 
ical states i4:i(1270) and i^i(1400). The rates of 
B — + KiAT^ decays, which are experimental inputs to 
the calculation of the bounds on |Aa|, must be derived 
from the measurement of the rates oi B ^ iiri(1270)7r 
and B — > Ki{1A00)tt decays. The BFs for these pro- 
cesses have recently been measured by BaBar f22|. 
The i4ri(1270) and i4ri(1400) axial vector mesons are 
broad resonances with nearly equal masses. In the fol- 
lowing, we will refer to them collectively as Ki. The 
ii'i(1270) and if i (1400) mesons decay to the same fi- 
nal state Ktttt, although through different interme- 
diate states. However, since the intermediate de- 
cays proceed almost at threshold, the available phase 
spaces overlap and interference effects can be sizeable. 
The analysis strategy relies on the reconstructed Kittt 
invariant mass spectrum in the [1.1, 1.8] GeV range to 
distinguish between ifi(1270) and ifi(1400), includ- 
ing interference effects in the signal model. 

A two-resonance, six-channel i^T-matrix model is 
used to describe the resonant Ktttt system for the sig- 
nal [1^. The production amplitude for channel i = 

{(K*Tr)s-.waveAK*7T)D-wave, pK , K^TT , foK , UjK} is 

given by 



= e'^^ ^(1 - »Kp)-.ip„ 



(23) 



where 6i are offset phases with respect to the {K*tt)s 
channel. 



faifaj fbifbj 



Ma - M Mb - M' 
and P is the production vector 



(24) 



p fpafai fpbfbi (0^\ 

Ma^M^ Mb-M- ^^^^ 

The labels a and b refer to ifi(1400) and ii:i(1270), 
respectively, and the indexes i and j refer to the final 



states of Ki decays. The decay constants fai, fbi, 
and the /^-matrix poles Ma and Mb are real. The 
elements of the diagonal phase space matrix p{M) 
for the process Ki 3-1-4, 3 5-1-6 have been 
approximated with the form 



25ij / 2m*m4 ^. 
P = -rf a/ ^ ^ ' (M - m* - m4 + »A), (26) 

where M is the mass of Ki, rui is the mass of 4, m* 
is the mean mass of 3 and A is the half width of 3. 

The parameters of K and the offset phases 8i are 
extracted from a fit to the data collected by the WA3 
experiment [23| for the intensity of the Kittt channels 
and the relative phases. For the fit to WA3 data a 
background term is included in the production vector. 

The decay constants for the uoK channel are fixed 
according to the quark model [23]. The production 
constants fpa and fpb are expressed in terms of the 
production parameters C = {'&t4>)'- fpa = cosz?, fpb = 
siniJe"^, where ■& G [0, 7r/2], e [0, 2tt]. 

Signal Monte Carlo (MC) samples are generated 
by weighting the {Kt:t:)t: population according to the 
amplitude where the term {KTnr\i) 

consists of a factor describing the angular distribution 
of the Kttt: system resulting from Ki decay, an ampli- 
tude for the resonant tttt and Kt: systems, and isospin 
factors. The BF of Ki loK is accounted for as a 
correction to the total selection efficiency. 

The BF and the production parameters d, (p for neu- 
tral and charged B meson decays to ifi(1270)7r + 
Xi(1400)7r are extracted via a ML fit to the kinematic 
observables rriESi AS, a J-"isher discriminant based on 
event-shape quantities, the if tttt invariant mass mKir-n 
and an angular variable. Background from B decays 
to iir*(1410)7r and non-resonant B decays to K*!^ 
and pKn are taken into account as separate compo- 
nents in the fit. The dependence of the signal prob- 
ability distribution in niKTr-n and selection efficiencies 
on the production parameters C, is described by means 
of non-parametric templates P(mif^^]i9, (/>). 

Each event is classified according to the invariant 
masses of the tt+tt^ and if+Tr" {K'^tt^) systems in the 
{Kf) decay for B^ {B~^) candidates: events which 
satisfy the requirement 0.846 < rriKTr < 0.946 GeV 
belong to class 1 ("if* band"); events not included in 
class 1 for which 0.500 < m^rTr < 0.800 GeV belong to 
class 2 ("p band"); all other events are rejected. 

For the i?° modes a likelihood scan is performed 
with respect to and (f). At each point, a simultaneous 
fit to the event classes r = 1, 2 is performed. Although 
for events in the "p band" the signal to background 
ratio is worse than that for events in the "if* band" , 
MC studies have shown that including those events 
in the fit helps in resolving the ambiguities in the de- 
termination of the parameter 0. For the B'^ modes, 
simulations show that, due to a less favourable signal 
to background ratio and increased background from 
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B decays, the analysis is not sensitive to (j). A value 
(j) — -K IS therefore assumed and the scan is performed 
only with respect to -d. At each point of the scan, a 
fit to "if* band" events only is performed. 

Figure [5] shows the distribution of Ai?, ruEs 
and mx-n-K for the signal events obtained by the 
background-subtraction technique sPlot [l^l . 

The experimental two-dimensional likelihood C for 
■d and (j) is convolved with a two-dimensional Gaus- 
sian that accounts for the systematic uncertainties. 
The resulting distributions in d and cf) are shown in 
Fig. [S] (the 68% and 90% probability regions are shown 
in dark and light shading respectively, and are de- 
fined as the regions which satisfy C{r) > Cmm and 
/£(r)>£„„„ ct))AMcj) = 68% (90%)). 

A combined signal for decays to ii'i(1270)+7r^ 
and ii'i(1400)"'"7r~ is observed with a significance 
of 7.5(7, while there's evidence for 5+ decays to 
is:i(1270)"7r+ and Ki{UmfTT+ at 3.2cr. The mea- 
sured BFs arc B{B^ X+tt" + K'+tt') = 3lt^ x 
10-6 and B{B+ -> K^t:+ + K'^t:+) = 29li? x 10"*^ 
(< 82x IQ-^ at 90% probability), including systematic 
uncertainties 22]. 

The probability distributions for the B 
Ki{127G)tt, B -> Xi(1400)7r, and B Kiatt 
BFs are derived by setting the production param- 
eters {fpa,fpb) equal to (0, e*"* sin?9), (cost?, 0), and 
(|/pa|cos6', -|/pA|sin6'), respectively, where fpA = 
cos cos O—e^"^ sin sin 9 and is the Ki mixing angle. 
A value 9 = 72° is used [H]. 

Including systematic uncertainties the following 
values are obtained (in units of IQ-^): B{B''^ 
is:i(1270)+7r-) = ITtli, >3{B" Xi(1400)+7r-) = 
171^, BiB° - K+A^-) = 14t?„, BiB+ ^ 
ii:i(1270)°7r+) < 40, B{B+ ifi(1400)°7r+) < 
39, B{B+ K%n+) < 36, where Jhe upper lim- 
its are evaluated at 90% probability [22j |. 

5.2. Extraction of a 

A MC technique is used to estimate a probability 
region for the bound on |Aq;|. The CP-averaged rates 
and CP-violation parameters participating in the es- 
timation of the bounds are generated according to the 
experimental distributions; a summary of the experi- 
mental values used as input to this calculation is pro- 
vided in Table [IVl 

For each set of generated values, the bound on |Aa| 
is evaluated. The limits on |Aa| are obtained by 
counting the fraction of bounds within a given value 
and the results are |Aa| < 11.1°(13.1°) at 68% (90%) 
probability (2^ . 

The angle a is extracted with an eight-fold ambi- 
guity in the range [0, 180]°. The eight solutions are 
a = (11±7±11)°, a = (41±7±11)°, a = (49±7±11)°, 
a = (79 ± 7 ± 11)°, a = (101 ± 7 ± 11)°, a = 
(131±7±11)°, a ^ (139±7±11)°, a = (169±7±11)°. 



Table IV Summary of the B ai(1260)7r and B 
ai(1260)-ftr branching fractions (in units of 10"®) and of 
the form factors (in MeV) used in the calculation of a. 



B(af7rT) [20] 


B{a-K+) [20] 


B(a+K°) [20] 


33.2 ±3.8 ±3.0 


16.3 ± 2.9 ± 2.3 


33.2 ±5.0 ±4.4 


U f4] 


U [4] 


/ai [25] 




130.4 ±0.2 


155.5 ±0.9 


203 ± 18 


207 ± 20 



Assuming that the strong phase S is negligible [2l|, 
only two solutions are still allowed. Considering only 
the solution consistent with the results of global CKM 
fits, a = (79± 7± 11)°. 

6. Conclusion 

Recent updates of measurements related to the de- 
termination of a have been presented. 

The first measurement of the branching fraction of 
B — !■ KiTT decays, combined with the input from the 
analysis of the time-dependent CP-violation asymme- 
tries in i?° — > ai(1260)^7r='= decays and of the branch- 
ing fractions of P ^ ai(1260)if decays, allows to 
measure a in the ai(1260)7r system. This novel de- 
termination of a is independent from, and consistent 
with, the current averages, which are based on the 
analysis of the tttt, pn, and pp systems only. 

With the new update of the B~^ p'^P'^ branching 
fraction and longitudinal polarization fraction mea- 
surements, the determination of a in the pp system 
has reached the unprecedented precision of 7%, com- 
parable with the 5.3% precision achieved in sin 2/3 
measurements . 

In the TTTT system, the updated measurement of CP- 
violating asymmetries in P" — s- tt+tt- decays provides 
a 6.7(7 evidence of CP violation. P" — > tt^tt" decays 
branching fraction and direct CP asymmetry are in- 
put to the isospin analysis of P — > tttt decays that is 
used to constrain the effect of penguin pollution on 
the extraction of a. 

All the measurements described in this work have 
been performed on the final BaBar sample. Most of 
them are still limited by statistics, and improvement 
may come from next generation, very high luminosity 
facilities. 
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Figure 8: sPlot projections of signal onto mss (left), AE (center), and rriKTrw (right) for B class 1 (top), B'^ class 2 
(middle), and B^ class 1 (bottom) events: the points show the sums of the signal weights obtained from on-resonance 
data. For niES and AE the solid line is the signal fit function. For mK-mr the solid line is the sum of the fit functions of 
the decay modes 7^1(1270)^ + Ki{lAm)-K (dashed), K*{1A1G)'k (dash-dotted), and if*(892)7r7r (dotted), and the points 
are obtained without using information about resonances in the fit, i.e., we use only the niES, AE, and variables. 




Figure 9: 68% (dark shaded zone) and 90% (light shaded zone) probability regions for i9 and (j) for the (a) B° and (b) 
_B+ modes. 
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